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Abstract. Functionals of Brownian bridge arise as limiting dis- 
tributions in nonparametric statistics. In this paper we will give 
a derivation of distributions of extrema of the Brownian bridge 
based on excursion theory for Brownian motion. Only the Poisson 
character of the excursion process will be used. Particular cases of 
calculations include the distributions of the Kolmogorov-Smirnov 
statistic, the Kuiper statistic, and the ratio of the maximum posi- 
tive ordinate to the minumum negative ordinate. 
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1. Introduction. 



Functionals of Brownian bridge arise naturally as as ymptotic distri- 

butions in nonparametric statistics. For an overview see lShorack and Wellner 
( 119861 ). The discussion in this paper will be limited to extrema of Brow- 
nian bridges even though the te chniques are app l icable to a wider class 
of homogeneous functiorials. Se e ICarmona et al.l (ll999[).|Pitman and Yor 
( 119991 ) , IPitman and Yorl (119981 ) and IPerman and Wellnerl (119961 ) . 

Let {Ut : < t < 1) be the standard Brownian bridge. Define 



(1.1) 

and 
(1.2) 



max Ut 

0<t<l 



— min Ut 

0<t<l 



m 



min{M+,M"} and M = max{M+, M"}. 



The distribution of M"*" was f irst computed by Smirnov (1939 ). For 

an elementary derivation see Govindarajulu and Klotj (Il973 ). The 

derivation of the distribution of M was given by iKolmogoro 3 (Il933h . 

Here we will give a derivation based on excursion theory for Brownian 
motion. It will be shown that the Laplace transform of M can be 
derived from that of M"*", and the distribution of can be derived 
by an elementary application of the reflection principle for Brownian 
motion. The Laplace transforms can be inverted explicitly as infinite 
series. In a similar fashion one can derive the joint distribution of 
(M+,M-). 

Let {Bt : if: > 0) be standard Brownian motion and denote the last 
exit time from before time t hy gt = sup{s < t : Bs = 0}. Further- 
more, define {L{t) : t > 0) to be the local time process at level for 
Brownian motion normalized so that 



:i.3) 



L(t) = max Bs 

0<s<t 



See iRevuz and Yorl (119991 ) for definitions 



The following lera ma is well known. See lLevyi (ll948l ). lDynkid (Il96ll ) 
Barlow et all (Il989f ). 



Lemma 1. The distribution of gi is Beta(l/2, 1/2). Given gi, the 
process {Bt : < t < gi) is a Brownian bridge of length gi, and 
the rescaled process {B{g\u) / ^fgi : Q < u < 1) is a Brownian bridge 
independent of gi. 

If Sq is an exponential random variable with parameter 9 indepen- 
dent of B define gsg = supju < Sg : Bu = 0}. Scaling properties of B 
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imply that the rescaled process 

B= (^^: 0<t<l^ 

is Brownian motion on the unit interval independent of Sg and the 
corresponding last exit time cji = sup{t < 1 : i?f = 0} is by Lemma 
[T]a Beta(l/2, 1/2) random variable independent of Sq. It follows that 
9Sg = 91- Sg is r(l/2, 9). Furthermore, given gsg, the process {Bt : < 
t < gse) is a Brownian bridge of length gsg and the rescaled process 
defined by 

U=(^^:0<t<l] 
V / 
is standard Brownian bridge on the unit interval independent of gsg 
and Sg. 

Hence if 7 ~ r(l/2, ^) is a random variable independent of the Brow- 
nian bridge U then 

(1.4) v^(M+,M-)4 max 5*,- min Bt) . 

\0<t<gsg 0<t<gsg J 

In Section 2 excursion theory will be used to derive the joint distribu- 
tion of the triple 

(1.5) ( max Bt,- min Bt,L{gsg)\ ■ 

\0<t<9Sg 0<t<gsg ) 

Only the fact that certain point processes are Poisson processes will 
be used rather than the more involved descriptions of the Ito excur- 
sion law. This will make it possible to derive distributions related to 
^/^(M"*", M~) by identifying this latter pair to the pair on the right in 
( 11.41) and conditioning on L (gsg)- In particular it will follow that given 
L (gsg) the random pair on the right in (11.41) are conditionally indepen- 
dent. Conditional independence will lead to formulae for the distribu- 
tion of the mimimum, the maximum, the sum, the difference and joint 
distribution of the random pair on the left in (11. 4p . The distributions 
related to (M+,M~) can be derived from those for ^ (M~^ , M~) by 
inverting Laplace transforms. 

2. Brownian excursions. 

Let B be standard Brownian motion. Let U = {w E C[0, 00), w(0) = 
0, 3R > 0, w{t) ^ iff t e (0, i?)} U {6} with 6 being the coffin state be 
the excursion space equipped with the cr-field generated by all coordi- 
nate maps. Let Tg = m{{u : L{u) > s} be the right continuous inverse 
of the local time process at level of i? and define the point process 



4 



MIHAEL PERMAN AND JON A. WELLNER 



e = ((s, Cs) : < s < L(S'e)) of excursions of the randomly rescaled 
Brownian bridge on the space (0, oo) x f/ by 

(2.1) e,(u) = llo<u<T,^Ts-)iu) Br^_+u 

for < u < Ts — when — > and = 5 if — Ts_ = where 
s < L{gs,). 

For w & U and w 6 define the duration of an excursion as 
R{w) = supj-u: w{u) 7^ 0} . 

Lemma 2. Let e be the point process defined in Ii2.1\) . 

(i) The random variable L{Se) is exponential with parameter \p2B . 

(ii) Conditionally on L{gsg) = t the point process e is a Poisson 
process in the space (0, t) xU with mean measure A x n ■ q-^^^'^) 
where A is the Lebesgue measure on (0,t) and n is Ito's excur- 
sion law for Brownian motion. 

(iii) Positive and negative excursions of {Bt: <t < gsg) are inde- 
pendent Poisson processes conditionally on L{gsg) = t. 

Proof: Let e be the excursion process of Brownian motion B. The 
excursion of Brownian motion straddhng Se can be interpreted as the 
first one in the process of marked Brownian excursions where an ex- 
cursion w is marked with probabihty 1 — e"^^*^""-* and marks are given 
independently. The marked and unmarked excursions are independent 
Poisson processes on (0, oo) x U with mean measures 

A X n ■ (1 - e-^^("')) and A x n ■ e^^^^"") 

respectively wh ere A is the Lebesgue rneasure and n is Ito's excursion 



law on U. See iRogers and Williamsl (Il987l ). p. 418 for details. The 



point process e is the initial portion of e up to but not including the 
first marked excursion. 

On the local time scale L{gsg) is the time of the occurrence of a 
point in a linear Poisson process hence exponential. By the reflection 
principle the density of maxo<s<( Bg is given by 

Hence the density of ma.xo<s<Se Bg is given by 

V2e 



IT 



e-^"/'H~'/^e-''dt = V2ee-'^'^ 







For a computation of this Laplace transform see lOberhettinger and Badiil 



( 119731 ) ■ , p. 41, formula 5.28. By (O]) this is also the density of L{Se). 
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The unmarked excursions up to L{gsg) are the excursions of the 
rescaled Brownian bridge {Bf-. < t < gsg)- It follows that given 
L{Sg) = t the process e is a Poisson process on (0,t) x U with mean 
measure A x n ■ e"^-^^'"^ where A is the Lebesgue measure on (0, t) and n 
is Ito's excursion law for standard Brownian motion. This proves (ii). 

Positive and negative excursions in the point process e are condi- 
tionally independent given L{Sg) = t hj independence properties of 
Poisson processes for disjoint sets which proves (iii). 



Let M"*", M~ and M be as defined in Section 1. Let U be Brownian 
bridge and 7 an independent r(l/2, 6) random variable. By the defini- 
tion of the Brownian bridge as Brownian motion conditioned to be at 
time t = 1 and an application of the reflection principle for Brownian 
motion one derives that for a; > 



P(M+ >x) = limP max Bt>x \ \Bi\ < e 

e^O \0<t<l ' 

= limP ( max Pt > xlPil < e I /P(|Pi| < e) 

e^O \0<t<l J I \< < 

= limP(2x-e < Pi < 2x + e)/P(|Pi| < e) 
(2.2) = 



It follows from the distribution of M"*" given by (12.21) that 

POO 

(2.3) P{y/^M+ >x) = y/ej^ / exp( -2x'^/s)s'^/^e~'^'ds 

= exp(-2xv^). 

(2.4) 



For a computation of this Laplace transform see lOberhettinger and Badii 
fll973[ ). , p. 41, formula 5.28. 

Turning to excursions recall that R{w) stands for the length of the 
excursion and denote = ma.Xuw{u). Define 



m(xi 



-eR{w) 



n{dw) . 



The point process of the maxima (e^: < s < L{Sg)) is a measurable 
map of the process e and hence conditionally on L{gsg) = t a Poisson 
process on (0, 00); see iKingman (Il993f ). The conditional mean measure 
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of the set [x, oo) is t ■ m{x). This imphes that 

(2.5) P ( max Bt > x\L{gs,) = t] = 1 - e'^''^^^) 

This in connection with Lemma [21 (i) and (12.31) gives that 
(2.6) 

poo 

Pi max Bt>x) = V2e / e"^*(l - e-*'"^^')) dt = 



'o<t<9Sg Jo y/2e + m{x) 

Comparing (12.31) and (12. 6p we obtain that 

(2.7) m(x) = ^ . 

Furthermore, given L{gsg) = t, by Lemma2, (iii), the positive and 
negative excursions of the point process e are conditionally independent 
Poisson processes. This implies for x,y > 
(2.8) 

P ( max Bt<x,- min Bt<y\ L{gsg) = t] = e~*'"(^)e-*'"(^) . 

\0<t<gSg 0<t<gsg J 

This also gives the joint distribution of the triple (11.51) . In particular 
(12.81) implies that two random variables 

(2.9) max Bi and — min Bt 

0<t<9Sg 0<t<gsg 

are conditionally independent given L^Se) = t. 

3. Examples of calculations 
3.1 Distributions of M and m. 

Let m and M be defined as in (11.21) . From (12. 6p we have 



V 2/9 

(3.1) P{^/lM <x) = = tanh(a;v/2^) 

\/26 + 2m(x) 



Let Fm be the cumulative distribution function of M. Writing out 
(13.11) . taking into account the independence of M and 7 and dividing 
both sides by VO we get 



(3.2) 



1 r°° 



tanh(xV26') 
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Oberhettinger and Badiil (jl973l ). p. 294, formula 8.51, give the inverse 



of the Laplace transform on the right of (13.21) as 



^^^^ V2: 



X 



oo 



5:(-irexp(-2nV/.). 
V2x vvrs ^ 



This yields 

oo 

Fm{z) = J2 (-I)"exp(-2n2z2), 

n=— oo 

or 

oo 

1 - Fm{z) = 2 exp( -2n'z'), 

n=l 

which is the formula for the distribution of the Kolmogorov-Smirnov 
test statistic. 

Turning to m = minlM"*", M~} observe that by (12. 8p 



poo 

(3.3) P i^m >x) = V2e / e-^* (l - e-*'"^"))' dt . 

Jo 

Integration yields 

(3.4) P{^m>x) = l- + ^ 



^29 + m(x) V29 + 2m{x) 

Denote by Fm the distribution function of m. From ( 13. 4p . independence 
of m and 7 and dividing both sides by ^/O we get 

1 , ^ , , M2 0.. tanh(a;\/2^) 26-^^^ 1 



The fi rst term has been inverted above, the second is given by lOberhettinger and Badii 
( Il973l ). p. 258, formula 5.87, and the third is elementary. Substituting 
z ioT x/ ^/s one gets 



1-Fm{z) = ^ (-l)"exp( -2712^2) + 2e" 

n=— 00 
00 

(3.5) = 2 J^(-l)"exp(-2nV). 



■2z 



2 



n=2 
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3.2 Joint distributions, sums, differences, quotients. 

Integrating out the conditioning variable in (12. 8p we get 
(3.6) P(77M+ < X, v^M" <y) = 

V2e 



29 + m(x) + m(y) 



^{x+y)V2e _ Q-{x+y)V2e 

sinh(xv^) smh{yV2e) 



smh{{x + y)V26) 
cosh((x + y)V2e) - cosh((x - y)V2e) 



smh.{{x + y)V29) 
coth((x + y)V2e) - cosh((x - y)V2e) csch((x + y)V2e) 



Let F{z,w) be the joint distribution function of the pair {M~^,M ). 

By (EH) 

(3.7) 

poo 

PiVlM+ < X, ^M- <y) = / F{x/^/~s, y/^/~s)s-'/^e'''ds . 

Jo 

The right side is given in (13.61) . Oberhettinger and Badii (1973), p. 
294, formula 8.52, give the inverse of the first term of the transform as 



1 



9. 



V2{x + y) V '2(x + i/) 



^ oo 

- E 



-2fc2(x+i/)2/s 



fc=— oo 

The inverse of the second term of the transform can be obtained from 



Oberhettinger and Badiil (Il973l ). formula 8.8.60: we find that the in- 
verse is 



1 



-9. 



(a;-l/)/2| 



V2ix + y) V x + y '2(a; + y)2 
-= 2^ exp -2{x + y) 



k=—oo 
oo 



+ * + i 

(x + y) 2, 



/s 



^ ^ exp(-2[A;(x + ?/) + x]Vs) . 



A;=— oo 
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Combining these yields 
P(M+ <z,M-< w) 



exp{-2k'^{z + wf) - ^ exp{-2[k{z + w) + z]'^) 



in agreement with IShorack and Wellnerl (119861 ). formula (2.2.22), page 
39. 

Another statistic whose distribution can be obtained from formula 
(12.81) is the Kuiper statistic defined as = M+ + M~ . It seems 
cumbersome to proceed from the joint distribution of M"*" and M~ so 
we use directly the formula (12.81) . Denote U = y^M+ and V = . 
The joint cumulative distribution function of V and V has been shown 
to be 



_ sinh(MV2^) sinh(t;V2^) 



smh{{u + v)V2e) 

A simple calculation shows that the cumulative distribution function 
oi U + V is given by 

PiU + V < z) = [ Guiu,z-u)du 
Jo 

where Gu is the partial derivative of G with respect to u. A straight- 
forward calculation yields 



, 2V29smh\(z-u)V29) 
Gu{u,z- u) = — 



sinh\zV2e) 
and integration gives 

z\/20 

P{U + V <z)= coth(2v^) - - 



sinh^(2v^) ' 

Using the fact that 7 and + M~ are independent we obtain the 
Laplace transform of the cumulative distribution function Fk of the 
Kuiper statistic as 



^ r F^(z/v/^)s-i/2e-^Ms = coth(zv^) 7^^^ . 

x/^io sinh2(zv^) 

After dividing by \/9 it remains to invert the two terms on the right 
and substitute for z/ y/s. The first term has been inverted above when 
deriving the joint distribution of M"*" and M~ . We get 
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1 1 °° 



k 

To invert the second term rewrite it as 



4^2. 



ze 



-22V26I 



?>mh^{z^) (1 - 



26^2 



A\f2z ^ ke- 



2kzV2e 



k=l 



for 2; > and ^ > 0. The inverse Laplace transforms of the terms in the 
sum are known. One can either recall that they are Laplace transforms 
of the hittin g time of level 2kz by standard BM, or use formula 5.85 
on p. 258 in lOberhettinger and Badiil (Il973l ) to show that for a > 



roc 

(3.8) / exp( -aV2s)e-^Ms = e"'^^ . 

Jo \'2ns^ 

Since all the terms are non-negative functions the order of summation 
and integration can be changed. Hence for ^ > the inverse Laplace 
transform of the second term is 



4^2^ y 4= exp( -2k''z^/s) = ^ V exp( -2k''z^/s) . 



k 

Substitute x = z/ ^/s to get 



(3.9) 

FAx) 



k=—oo 



-2k^x 



2^2 



- E « 

k=—oo 



2^2g-2/c2x^ 



2fc2x2 



fc=— oo 



Remark: IVervaatI (Il979l ) gives a construction of standard Brownian 
excursion from a Brownian bridge. Let U he a Brownian bridge [0, 1] 
and let a be the time when U attains its minimum on [0, 1] [a is a.s. 
unique). Then the process (e(t) : < t < 1) defined by 



e(t) = Ufj+t {mod 1) — U^it) 

is a standard Brownian excursion. It is a simple consequence of this 
transformation that the Kuiper statistic has the distribution of the 
maximum of the standard Brownian excursion and 03.91) is another 
derivation of the distribution of this maximum. 

For the difference U — V a straightforward computation yields for 

z>0 
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(3.10) P{U-V>z) = j Gu{u,u-z)du 

provided P{U > 0) = 1 and P(y > 0) = 1 which is the case for the 
variables in question. Using the joint cumulative distribution function 
yields 

sinh2(v^(M - z)) 



Gu(u,u - z) = 2V2e ^ ^ 
^ ^ sinh2(v^(2M-z)) 

FromdSinD it follows 

sinh2(v^(M - z)) 
sinh2(v^(2M-z)) 



p{u-v>z) = 2V2e I , ^ / ^d« 



, ^ , sinh2(v^M) 
2V2e / — ^ du. 



smh\V29{2u + z)) 

Since 7 and {M~^,M~) are independent the above identity, after can- 
celling \/9 can be rewritten as 



^P(M+-M- > ^)e^^'ds = 2 ^ A, ^ du . 

\^ Jo Vs Jo sinh^(v^(2u + 2)) 

This gives the Laplace transform of P{M~^ — M" > zj^^j^ as a 
function of s for fixed z. 

To invert this Laplace transform note that for fixed m > and 2; > 
the integrand on the right as a function of ^ is a Laplace transform. 
The inversion proceeds by first inverting the integrand for fixed u and 
then integrating with respect to u. For 6' > and -u > write 

sinh2(v^M) g2v^« _ 2 + e-2v^" 



sinh2(v^(2u + z)) e2v^(2"+-) - 2 + e-2^(2«+^) 

/'g2V29« — 2 + e~2v^«)A . g-2v^(2M+2) 



\ _ 2eV26'(2u+^) _|_ g-4V2e(2M+2) 

00 

^^2V2eu _ 2 + g-2V2e«^ ^ ^ g-2A:V2e(2u+^) 



fe=l 

The series expansion follows from noting that for < 1 

00 

^ ^ fc=l 
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and setting x = e 2v2e(2n+2)_ Multiplying out one obtains 
(3.11) 

oo oo oo 

^g-(2fc2+(4fc-2)u)v^ _ 2 ^ j^^-2k{2u+z)V2e _^ ^ ^g-{2fc2+{4fc+2)«)v^ 

fc=l fc=l A:=l 

All the terms in the sums can be inverted using (13.81) . From left to 
right in (13. lip the inverses of the sums are 

(3.12) -jL= V k{2kz + (4m - 2)) ■ e-(2'=-+(4'=-2)-)V2. ^ 

V27rs2 



fe=i 



^ oo 

(3.13) - 2-== Y.'^k\2u + z) ■ e 



{2k{2u+z)Y /2s 

3 ' 
k=l 



and 



(3.14) ^1= ^ ki2kz + (4n + 2)) ■ e-(2^-+(4^'+2)«)V2s _ 

The sums are absolutely integrable so the order of summation and 
inversion can be interchanged. To derive the distribution of M"*" — M~ 
multiply (I3.12p . (I3.13P and (13.140 by 2v^27r and integrate with respect to 
u. Again the change in order of summation and integration is justified 
by absolute integrability and Fubini's theorem. Summing up yields 
(3.15) 



-. / oo „ , oo oo 

_L I \^ -2^2^2/5 _ Sr^ -2k'^z^/s , \ ^ 



2^ ^--2k'^z'^/s 



^/il^4A;-2 ^ ^ Ak + 2 

^ \k=l k=l k=l 

Finally note that (13.15^ gives 

1 1 °° 1 



or 

oo ^ 

P{M^-M->z) = J2^^^-e 

k=l 

For a different approach see Kosorok and LirJ ( 1999 ). Turning to quo- 
tients define Q = /M~ . We can multiply the numerator and de- 
nominator by and choose 9 = 1/2. Conditionally on //(Se) = t an 
elementary calculation gives the conditional distribution function of Q 

as 

POO 

FQ\LiSe)=t{z) = - / tm'(x)e-*™(^')e-*"^(^^) dx . 
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Unconditioning and changing the order of integration gives 

m'ix) dx 



(3.17) 



Fa z) 



lo (1 + m(x) + m(xz))^ 
Substituting (12. 7p and observing that 

1 



m'(x) 



2sinh X 



gives 



Fo(z) 



1 /-^ (1 - e-2-)2(l - e-2-^)2 

2 sinh2x(l -e-2^(^+i)) 
sinh^zx dx 



dx 



sinh {{z + l)x) 



z + l 



1 - 



TTZCtg (^) 
Z + l 



where the last integral is given in iGradshteyn and RyzhikI (Il994l ). for- 
mula 3.511.9. 



3.3 COVARIANCE AND CORRELATION. 

Yet another functional of interest is the correlation between and 
M~ . By independence of 7 and (M^,M~) it will suffice to compute 
the covariance between U = and V = . Starting from 

(3.18) 

cov(f/, V) = E (cov(f/, V\L{gs,)) + cov {E{JJ\L{gs,)). E{y\L{gs,)) . 

we first notice that by (II. 4p . (II. 5p and (12.80 the variables U and V 
can be taken to be conditionally independent, given L{gsg)- The con- 
ditional covariance in (13.180 is 0. By independence and symmetry 
E{U\L{gs,)) = E{V\L{gs,)). One gets 

cov(f/,\/)=var(E(f/|L(^75j) • 
The conditional expectation is computed from (12. 8p as 



poo 

(3.19) E{U\L{gs,) = t) = (1 - e-*™(^)) dx 

^0 
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Denote ip{u) = E{U\L{gsg) = u). Note that for 29 = 1 the local time 
L{gsg) is exponential with parameter 1 and compute 

Ei^'iLigsg))) = 



POO 

dx / dy 



1 



1 



OO poo 



dx / dy 



1 + m{x) 1 + m(y) 1 + m(x) + m{y) 



-2t rt 



dt- 



1 - e 

OO g^2t 

1 - e-2* 



-2t 



1 _ e-2(^+j/) 



(2t-(l-e-'*))dt 



/■OO 1 

-^0 k=l 



/£ = ! 

7r2 1 
12 ~ 2 



The second line follows from Fubini's theorem and the third by chang- 
ing variables to x + y = t,y = v. From (12.81) we have E{U) = 
E{i){L{ge)) = 1/2 and hence 

(3.20) coy{U,V) = Y^T{^{L{gsg))) 

7r2_3 
~ 12 4 

From one derives E{M+) = E{M-) = y27r/4 and var(M+) = 
1/2 — 7r/8so from 

cov(?7, V) = E{-f) E{M+M-) - E'^{^) E{M+)E{M-) 

we have 

71^ 1 

E{M+M-) = - 2 • 

and 

(3.21) cov(M+,M-) = ^-i-^ 
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yielding 

TT^ 1 TT 

(3.22) corr(M+, M~) = \ ^ = -0.654534 . 

2 ~ 8 



16 



MIHAEL PERMAN AND JON A. WELLNER 



References 

Barlow, M., Pitman, J. and Yor, M. (1989). line extension mul- 
tidimensionnelle de la loi de Tare sinus. In Seminaire de Probabilites, 
XXIII, vol. 1372 of Lecture Notes m Math. Springer, Berlin, 294-314. 

BlANE, P. and YoR, M. (1988). Sur la loi des temps locaux browniens 
pris en un temps exponentiel. In Seminaire de Probabilites, XXII, 
vol. 1321 of Lecture Notes in Math. Springer, Berlin, 454-466. 

Carmona, p.. Petit, F., Pitman, J. and Yor, M. (1999). On the 
laws of homogeneous functionals of the brownian bridge. Studia Sci. 
Math. Hungar. 35 445-455. 

Dynkin, E. B. (1961). Some limit theorems for sums of independent 
random variables with infinite mathematical expectations. In Select. 
Transl. Math. Statist, and Probability, Vol. L Inst. Math. Statist, 
and Amer. Math. Soc, Providence, R.I., 171-189. 

Govindarajulu, Z. and Klotz, J. H. (1973). A note on the asymp- 
totic distribution of the one-sample Kolmogorov-Smirnov statistic. 
Amer. Statist. 27 164-165. 

Gradshteyn, I. S. and Ryzhik, I. M. (1994). Table of integrals, 
series, and products. Russian ed. Academic Press Inc., Boston, MA. 
Translation edited and with a preface by Alan Jeffrey. 

Kingman, J. F. C. (1993). Poisson processes, vol. 3. The Clarendon 
Press Oxford University Press, New York. Oxford Science Publica- 
tions. 

KoLMOGOROV, A. N. (1933). Sulla determinazione empirica delle 
leggi di probabilita. Giorn. 1st. Ital. Attuari 4 83-91. 

KOSOROK, M. R. and LiN, C.-Y. (1999). The versatility of function- 
indexed weight-rank statistics. J. Amer. Statist. Asso. 94 320-332. 

Levy, P. (1948). Processus Stochastiques et Mouvement Brownien. 
Suivi d'une note de M. Loeve. Gauthier-Villars, Paris. 

Oberhettinger, F. and Badii, L. (1973). Tables of Laplace trans- 
forms. Springer- Verlag, New York. 

Perman, M. and Wellner, J. A. (1996). On the distribution of 
Brownian areas. Ann. Appl. Probab. 6 1091-1111. 

Pitman, J. and Yor, M. (1998). Random brownian scaling identities 
and splicing of bessel processes. Ann. Probab. 26 1683-1702. 
URL [http :// stat . berkeley . edu/ user s/pitman/590 .pdf | 

Pitman, J. and Yor, M. (1999). Path decompositions of a brownian 
bridge related to the ratio of its maximum and amplitude. Studia 
Sci. Math. Hungar. 35 457-474. 



AN EXCURSION APPROACH TO MAXIMA OF THE BROWNIAN BRIDGE 17 

Revuz, D. and YoR, M. (1999). Continuous martingales and Brow- 
nian motion, vol. 293 of Grundlehren der Mathematischen Wis- 
senschaften [Fundamental Principles of Mathematical Sciences]. 3rd 
ed. Springer- Verlag, Berlin. 

Rogers, L. C. G. and Williams, D. (1987). Diffusions, Markov 
processes, and martingales. Vol. 2. Wiley Series in Probability and 
Mathematical Statistics: Probability and Mathematical Statistics, 
John Wiley & Sons Inc., New York. Ito calculus. 

Shorack, G. R. and Wellner, J. A. (1986). Empirical processes 
with applications to statistics. Wiley Series in Probability and Math- 
ematical Statistics: Probability and Mathematical Statistics, John 
Wiley & Sons Inc., New York. 

Smirnov, N. (1948). Table for estimating the goodness of fit of em- 
pirical distributions. Ann. Math. Statistics 19 279-281. 

Smirnov, N. V. (1944). Approximate laws of distribution of random 
variables from empirical data. Uspehi Matem. Nauk 10 179-206. 

Smirnov, N. V. (1947). Sur un critere de symetrie de la loi de dis- 
tribution d'une variable aleatoire. C. R. (Doklady) Acad. Sci. URSS 
(N.S.) 56 11-14. 

Vervaat, W. (1979). A relation between Brownian bridge and Brow- 
nian excursion. Ann. Probab. 7 143-149. 

MiHAEL Perman, Institute for Mathematics, Physics and Mechanics, 
Jadranska 19, SI-1000 Ljubljana, Slovenia. 
E-mail address: mihael.permanSuni-lj .si 

Jon a. Wellner, Department of Statistics, University of Washing- 
ton, B320 Pedelford Hall, Seattle, WA 98195-4322, U.S.A. 
E-mail address: jaw@stat.washington.edu 



